arXiv:1505.02266v2 [cond-mat.soft] 27 Jul 2015 


cpi uiaii 


A modified Poisson-Boltzmann theory: Effects of co-solvent po¬ 
larizability 

Yu. A. Budkov 1 ’ 2 ( a ), A. L. Kolesnikov 3,4 and M. G. Kiselev 1 

1 G.A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, Ivanovo, Russia 

2 National Research University Higher School of Economics, Department of Applied Mathematics, Moscow, Russia 

3 Institut fur Nichtklassische Chemie e.V., Universitat Leipzig - Leipzig, Germany 

4 Ivanovo State University, Ivanovo, Russia 


PACS 82.45. Gj - Electrolytes 

PACS 33.15.Kr - Polarizability of molecules 

PACS 65.40. gk - Electrochemical properties 

Abstract -In this paper within a field-theoretical approach taking into account explicitly a co¬ 
solvent with a nonzero dipole and a polarizability tensor, we derive a modified Poisson-Boltzmann 
equation. Applying the modified Poisson-Boltzmann equation, we formulate a generalized Gouy- 
Chapman theory for the case when an electrolyte solution is mixed with a polar co-solvent having 
a large polarizability. We show that an increase of the co-solvent concentration as well as the 
co-solvent polarizability lead to a significant increase of differential capacitance at sufficiently 
high surface potentials of the electrode, whereas the profile of the electrostatic potential becomes 
considerably more long-ranged. On the contrary, an increase in the permanent dipole of the 
co-solvent only weakly affects the differential capacitance. 


Introduction. — Up to now the Poisson-Boltzmann 
(PB) equation remains the simplest and the most effective 
theoretical tool for describing distribution of charged par¬ 
ticles near charged surfaces in biophysics, colloid chem¬ 
istry, electrochemisrty, etc. However, the PB equation 
contains many assumptions that make it inapplicable to 
real physical systems. The mean-field nature of the PB 
equation does not allow us to take into account the ef¬ 
fects of the ionic correlations, whereas the consideration 
of the solvent as a continuous dielectric medium makes 
it impossible to study the effects of the solvent molecular 
structure. These two factors motivated the researchers to 
find new approaches the PB equation modifications during 
the last two decades 00 - 

Nowadays, great efforts have been made to modify the 
PB equation with respect to ionic correlations [I)|8], the 

11:16] , polarizability and 


dipole structure of the solvent 9 


permanent dipole of ions 13 -16 as well as their excluded 
volume [17)421], and finally solvent quadrupolarizability 
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chemical characteristics of real solutions, such as dielectric 
permittivity and activity coefficients 11 22 . 

Nevertheless, despite the evident success in extending 
the pure PB-paradigm there is still no answer to the ques¬ 
tion whether or not it is possible to incorporate into the 
PB equation the effects of the solvent polarizability in the 
electric field of the ions and external charges. 

Schroder and Steinhause presented 


25 a research of 


an influence of cation polarizability on the structure and 
dynamics of an ionic liquid within Molecular Dynamics 
(MD) simulations. It is shown that when the cation po¬ 
larizability increases, the structure of ionic liquid remains 
almost unchanged. On the contrary, the rotational and 
translational diffusion sufficiently increase. Recently, us¬ 
ing MD simulations for the set of ionic liquids was shown 
that cation polarizability has a significant effect on a shape 
of the lines of a vibration density of states 26 . These ob- 


22 . These modifications sufficiently improved the theory 


of electrolytes that enabled calculations of the physico- 


servations indicate that polarizability significantly affects 
the characteristics that are mainly determined through the 
long-range correlations of the solution particles. The dif¬ 
ferential capacitance of an electric double layer is one of 
the physical values that is mainly determined by the long- 
range correlations. Thus, to correctly evaluate the differ¬ 
ential capacitance and the electrostatic potential at long 
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distances from the charged electrode, it should be impor¬ 
tant to incorporate the polarizability effect of the solvent 
and the electrolyte ions. Moreover, simultaneous incorpo¬ 
ration of the polarizability effects of the solvent and ions 
into the PB equation is especially important for electro¬ 
chemical applications, where one of the main problems is 
the search factors that can significantly increase a magni¬ 
tude of the differential capacitance. 

In particular, two questions should be answered: 

• How much the differential capacitance will change 
when the polarizability of the solvent will be taken 
into account? 


• What is the difference between the effects of inclusion 
a permanent dipole and a polarizability of the solvent 
on the differential capacitance? 


Addressing these two questions, we derive within the 
field-theoretical approach a modified PB equation taking 
into account explicitly a solvent with a nonzero permanent 
dipole and a polarizability tensor. Using this modified 
PB equation, we formulate a generalized Gouy-Chapman 
theory 23 24 for the case when the electrolyte solution 
is mixed with a polar co-solvent having a strong polariz¬ 
ability. We calculate the differential capacitance and the 
electrostatic potential profile near the charged electrode at 
different co-solvent concentrations. We show that adding 
a strong polarizable co-solvent to the electrolyte solution 
leads to a significant increase of the differential capaci¬ 
tance under the sufficiently large surface potentials of the 
electrode and to a more long-ranged electrostatic poten¬ 
tial profile. We obtain that with the increase of the co¬ 
solvent polarizability, the differential capacitance dramat¬ 
ically increases. On the contrary, the permanent dipole of 
the co-solvent has only a weak effect on the value of the 
differential capacitance. 


Theory. — We consider an electrolyte solution con¬ 
taining N. |_ point ions carrying a charge z + e ( z + > 0), N_ 
ions carrying a charge Z-e (z- < 0 ), and a solvent which 
we shall model as a continuous dielectric medium with di¬ 
electric permittivity e s . In addition, we also consider N c 
molecules of a co-solvent which has a permanent dipole p 
and a polarizability tensor a. Moreover, we assume that 
the positively charged ions (cations) have also a polariz¬ 
ability tensor d + . The latter assumption covers the case 
when the electrolyte is an ionic liquid. The discussed sys¬ 
tem confined to a volume V and has the temperature T. 
Since in this work we shall discuss the effects of polarizabil¬ 
ity and permanent dipole only, for the sake of simplicity 
we neglect the excluded volume of the ions as well as the 
co-solvent molecules. The variants of models in which the 
excluded volume of the ions has been taken into account 


have been developed in recent works 17 21 


A canonical partition function of the above-mentioned 
system can be written in the following form 


Z= dT c 


dTie 


~(3H 


( 1 ) 


where the following notations for the integrating measures 
have been introduced 


' , A- 3N +A~ 3N - ' N - 

dTi(..)= + 


. N- 


N+\N-\ 


n dr t dv k / n dr i (••) 


k =1 


1 =1 


( 2 ) 


and 


A, 


-3N C r N ° 


N 0 


dr c (..) = I n d PkdV k I n **(••)■ ( 3 ) 


k =1 


k= 1 


The Hamiltionian of the system can be written as 

1 Nc 1 N+ 

H = 2 E^^ + 2 yp, + a^, + 

k =1 1=1 

+ \j dr J dr'p(r)V c (r- r')p(r'), (4) 

v v 


where V k and are fluctuating dipoles of k- th molecule 
of the co-solvent and Z-th cation, respectively; a and d + 
are the polarizability tensors of the co-solvent molecules 
and the cations, respectively. Moreover, a local charge 
density p(r) of the system has been introduced: 


N + N- 

p( r) = z+ e J2H r ~ r i)+ z - e J2 s ( r ~ r k) 

k =1 k= 1 

N c N + 

-XE Pfe Vl5 ( r ~ r k)~^2' P k VS ( r ~ r fe)+ZW( r), 

k =1 k =1 

(5) 

where pfc is the permanent dipole of the fc-th co-solvent 
molecule; p e xt{r) is the density of an external charge; 
V c (r — r') = g 1^71 is the Coulomb potential. It should 
be noted that in Q for the sake of simplicity the electro¬ 
static self-energy has been omitted. 

Using the standard Hubbard-Stratonovich transforma¬ 
tion, after a calculation of the Gaussian integrals over the 
variables Vj and and an integration over the orienta¬ 
tions of the dipoles pfc, we arrive at 


f Vlfi f dr(Vip(r)) 2 +if drp e:rt (r)ip(r) 

Z = Z id>i Z id , c / -f e 

J ^0 

( 6 ) 




N_ 


, r )) 2 

where C 0 = J Vipe v 
constant; /? = 1 /fcsT; 


is the normalization 


A ~3N + . -3iV_ yN + +N_ 

Z ld ,i = - (2nk B Tdeta+) N+/2 (7) 


is the ideal partition function of the ions, 


^-3N c yN a 

Zid, c = -(2 nksT det d) ^ /2 ( 8 ) 


N r \ 
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is the ideal partition function of the co-solvent molecules; 
A± and A c are the thermal de Broglie wavelengths of the 
ions and the co-solvent molecules, respectively. The single¬ 
particle partition functions of the ions and the co-solvent 
molecules that are immersed into an auxiliary fluctuating 
electric field with potential iip( r) take the following forms: 


_ I ^giz+ev(r)-kf-fVv>(r)a + V^(r)^ 


Q-~[ 


/*(■ 


n N- 


,iz-e<p( r) 


1 + 7 -' e 

v 

P- f dr(e* z ~ e ‘ p ^ — l) 


- 1 


M P+ f dr ( e *~ + 

Q+ + [<p]*e - 


ZV>( r )~ V^( r )«+ V^( r )_i 


Q 


N c 


J>*S dr^e-^ - 1 ) 


^ p.cr. — 


Z = f 

Zid,iZid,c J ^0 


Then, we calculate the functional integral (15) within a 
standard saddle-point approximation [2j. An equation for 
the saddle-point has the form 


M r ) 


= 0 . 


(18) 


(9) 

Q4<P] = Jye i *~ e *U, ( 10 ) 

V 

QM = / | e -^v, Ma v, H smp|Vy(r)| ^ 

J v p|V<p(r)| 

v 

In the thermodynamic limit (N± —> oo), we have the 
following asymptotic relations 27 


Calculating the variational derivative in (18) and intro¬ 


ducing a notation for an electrostatic potential ip(r) = 
—iksTiffc), one can obtain the following equation 


z_etp(r) 


V (e(r)VV’(r)) = — 4:TTp ext (r) — Anz-ep-e 


— Anz+ep+e k B 1 


z + eip( r) Vt/)(r)Q:_|_ Vi/>(r) 


(19) 


where a tensor of the local dielectric permittivity has been 
introduced 

vy-(r)avy.p) sinh ftp\\7ip(r)\ 
e(r) = £ S I + 47 rp c e 2k s T 


PP I V-0(r)| 


( 12 ) 


(13) 


a 


P 2 L(f3p\\7ip(r)\) 


k B T /3p|VV’(r)| 


+ 4irp+a+e 


z + eV»(r) [ V^(r)a + Vi/i(r) 


( 20 ) 


where p± = N±/V are average concentrations of the ions 
in the bulk solution. 

Analogously, in the thermodynamic limit {N c — > oo) we 
can obtain 


where L(x ) = coths — 1/x is a Langevin function. The 
tensor of dielectric permittivity in the bulk solution (where 
ip = 0 ) is determined by the following expression 


£ = £ S I + Anp c d + I 


3 k B T 


+ 4tt p+a +1 


( 21 ) 


(14) 


Thus, the excess partition function of the system takes the 
form 


(15) 


where the following short-hand notation for the integrand 
has been introduced 

<%] = ^ J dr (V(p(r )) 2 - i J drp ext (r)<p(r) - W^cp], 
v v 

(16) 

where a ’’functional of interaction” Wj[p\ has been also 
introduced 

Wj[p\ = p c [ dr f e ~¥^(»v y (r) sinp|Vy(r)| _ A 

J \ P |Vp(r)| J 

v 

+ p + J dr ^e iz+ev( - r )~^~ v ^( r )“+ v vd r ) — i^ 
v 

+ P- J dr (e iz - ev(r) - l) . (17) 


where I is the identity tensor; p c = N c /V is a bulk co¬ 
solvent concentration. 

The equations ( I9p0 ) determine the electrostatic po¬ 
tential ip{ r) and the tensor of local dielectric permittivity 
within a self-consistent field approximation. The latter is 
a generalization of the classical PB theory for the case of 
an explicit account of the polar polarizable co-solvent. If 
one assumes that d = d+ = 0 and e s = 1 , then one arrives 


at the expression which is obtained in the work 111 


e ( r ) = 1 + G(/3p\Vil’(r)\), 


( 22 ) 


where G{x) = s -^L(x). 

In the case of nonpolar co-solvent (p = 0), we obtain 
the following expression for the tensor of local dielectric 
permittivity: 


e(r) = £ S I + Anp c 
+ 47rp + d+e 


Vt/>( r)aVi/i(r) 

de 2 k B T 

z_|_e^(r) | Vt/>( r)a + 7i/>(r) 


(23) 


Finally, an excess electrostatic free energy of the solu- 
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tion can be written in the following form 

FelbP] =-^ J dr(Vip(r)) 2 


eV'(r) 


P-k B T J dr |^e k B T —1 
v 

/• f ~ + eiHr) Vy J (r)a + V^(r) 

p + k B T / dr ( e k B T + 2k B T — 1 


p c k B TJdrye (3p\Vip(r)\ 1 

V 


drp ext (r)ip(r). 


(24) 


The differential capacitance of the electric dou¬ 
ble layer: Effect of co-solvent polarizability and 
permanent dipole. — As an illustration of the applica- 

we formulate the 


201 


23 


24 . We consider 


— (e(z)ip'(z)) = 8npesmh/3eip(z), 
dz 


(25) 


e(z) = e s + 47 rp c e 


M'(z) 


24 28 


— e(O) , 0 / (O) = 47rcr. 


(27) 


tion (251 


e s £ 2 (z) 

87T 


£{z) 


J D(£')d£' = 2 pk B T (cosh(/3 eip(z)) — 1 ), 

(28) 


where £(z) = — ip'(z ) is the electric field and the auxiliary 
function D (£) takes the form 


see eq. (29) below 


In the case of a nonpolar solvent (p = 0), the first inte¬ 


gral (28) can be expressed through the elementary func¬ 


tions as 

e s £ 2 {z) 

8n 


+ 


*£*( = ) \ °£ 2 (.Q 

p c k B T 1 — e 2k B T ) + p c a£ 2 (z)e 2k s T 


= 2pk B T (cosh fleip(z) — 1). 


(30) 


To integrate the equation (28 1 , we should first solve it as 


tion of the modified PB equation (19 
generalized Gouy-Chapman theory 
a system containing a charged electrode, which we shall 
model as a charged flat surface with a fixed surface charge 
density a, the point ions of 1:1 electrolyte (z + = —Z- = 1 ), 
and the point molecules of the polar polarizable co-solvent 
with a polarizability a = al and a permanent dipole mo¬ 
ment p. In this case the average concentrations of ions in 
the bulk are equal, i.e. p + = p- = p. For the sake of 
simplicity, we assume that = 0 . 

Choosing z axis perpendicular to the electrode and plac¬ 
ing the origin on it, one can write the modified PB equa¬ 
tion in the following form 


where the effective dielectric permittivity e(z) takes the 
form 

**'(*) 2 sinh f3pil>'(z) 


a transcendental equation with respect to £ = —ip'(z ) at 
different values ip. Thus, we obtain a function £ = £(ip). 
The next step consists of solving the equation ip' = —£(ip). 

Before calculating the electrostatic potential profile 
ip(z), we estimate the differential capacitance C. As it was 
already pointed out in Introduction, this quantity should 
be sensitive to an account of the co-solvent polarizabil¬ 
ity. The differential capacitance can be determined by the 
following relation 

C=4^, (31) 

dipo 

where ipo = ip( 0 ) is the surface potential of the electrode, 
which in contrast to the surface charge density o is usually 
an experimentally controllable parameter. 

To perform the numerical calculations, we introduce di¬ 
mensionless parameters in the following way 

P = p/el B , a = a/l 3 B £ s , 0 = l 3 B p, 7 = l B p c , (32) 

where l B = e 2 /£ s k B T is the Bjerrum length. More¬ 
over, we introduce the dimensionless electrostatic poten¬ 
tial u = /3eip, the electric field £ = £/3el B , and the dis¬ 
tance from the electrode z = z/l B . We also introduce 
the dimensionless surface charge density d = crftel B /£ s , 
which is related to the dimensionless electric field at the 


electrode £q = £(0) via the boundary condition (27). The 


( p 2 L(j3pip'(z))\ 

k ( a + kPf^4W)- (26) 

In order to obtain the electrostatic potential profile ip(z), 
we have to formulate the standard boundary condition 


latter can be rewritten in the dimensionless form as 


1 + 47T7e 2 £(£q) ) £0 = 47tct, 


(33) 


where £(£ 0 ) = (d+pg^). 

Further, using the definition of the reduced differential 
capacitance 

(34) 

CLUq 


Using the standard method of integrating the equations 
of this form, taking into account the second boundary con¬ 
dition ip'( 00 ) = 0 , we can obtain the first integral of equa- 


equation (28) at z = 0 is written in the dimensionless form 

£0 

£q + 8n J dvD (v) = I 67 t6 (coshito — 1 ), (35) 

0 

after some algebra, one can obtain 

26 sinh ug 


C = 


£o(uo) 


(36) 
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D{£) = p c e 3 


. sinh (3p£ 
Pp 



p 2 L(/3p£) \ 
k B T j3p£ ) 


a£ 2 

ksT 


+ (3p£ coth f3p£ 


P 

k B T 


L'{Pp£) 


L(PpS) \\ 

Pp£ ))■ 


(29) 





Fig. 1: (a) The differential capacitance C as a function of the 
surface potential uo at different 7 and at fixed parameters a = 
10 —3 , 6 = 10 -3 and p = 0. (b) The differential capacitance as 
a function of the surface potential u 0 at the different co-solvent 
polarizability a and at fixed 7 = 0.1, 9 = 10 -3 and p = 0. 


where the function £q = £o(uo ) is determined implicitly 
through ( [35| ) and Uq = w(0). The dimensionless auxiliary 
function D (v) has the following form 


see eq. (37) below 


At 7 —> 0 we obtain £q ~ 4-\/27r 9 sinh ^. In this case, 
we arrive at the well-known result obtained by Grahame 
29 in the framework of the standard Gouy-Cliapman 


model 


C = 



(38) 


Fig la shows the differential capacitance as a function of 
the surface potential at different values of parameter 7 and 
fixed values a = 10 -3 , 9 = ICG 3 , and p = 0 which approx¬ 
imately correspond to the real systems (such as benzene or 
toluene, dissolved in a dilute aqueous electrolyte solution). 
As it can be seen from Fig. 1, the differential capacitance 
does not almost depend on the co-solvent concentration 
at the small values of uq. However, adding the polarizable 
co-solvent at sufficiently large surface potential leads to 


Fig. 2: The dependences of the differential capacitance C on 
the surface potential uo at different permanent dipole p of the 
co-solvent. The data are shown for 7 = 0.1, 6 = 1CP 3 , a = 0 . 


a significant increase in the differential capacitance. The 
same trend takes place at increasing the dimensionless co¬ 
solvent polarizability (see, Fig. lb). Fig.2 demonstrates 
the dependences of the differential capacitance on the sur¬ 
face potential at different values of the dimensionless per¬ 
manent dipole p of the co-solvent molecules. In contrast 
to the previous case, when the polarizability increased, 
an increase of the co-solvent permanent dipole leads only 
to a slight increase of the differential capacitance. The 
dependences of the electrostatic potential profiles u(z) at 
different 7 are marked in Fig. 4. As it can be seen, the 
presence of a polarizable co-solvent in the electric double 
layer makes the electrostatic potential considerably more 
long-ranged. 


Summary. — In this letter, in the framework of a 
field-theoretical formalism we have derived a modified 
Poisson-Boltzmann equation for the case of an explicit 
account of a polarizable polar co-solvent. Moreover, an 
expression for the tensor of the local dielectric permit¬ 
tivity has been obtained. As an important application 
of the developed mean-held theory, we have formulated 
a generalized Gouy-Chapman model of the flat electric 
double layer with the explicit account of the added po¬ 
lar polarizable co-solvent. We have shown that adding 
the strong polarizable co-solvent leads to a significant in¬ 
crease of the differential capacitance in a region of the 
sufficiently high surface potentials of the electrode. We 
have also shown that increasing the co-solvent polarizabil¬ 
ity leads to a significant increase in the differential capac¬ 
itance, while growth of the permanent dipole gives only a 
slight enhancement the latter. As it follows from (35-36) 
the differential capacitance is proportional to the value of 
the local dielectric permittivity at the electrode. Thus, it 
should be sensitive to a change of the co-solvent concentra¬ 
tion as well as the co-solvent polarizability (see, (26)). It 
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~ , . sinh pv 

D(v) = ye 2 --— 


at + p 


L(pv) 


(av 2 + pv coihpv) + p 2 I L'(pv) — 


Ljpv) 

pv 


(37) 



Fig. 3: The electrostatic potential profiles u(z) at different 
values of the parameter 7 . With the increase of co-solvent 
concentration the electrostatic potential becomes substantially 
more long-ranged. The data are shown for a = 0.005, 9 = 
10 _i , Mo = 10 , p = 0 . 


explains such dramatic behavior of the differential capaci¬ 
tance when these variables increase. In addition, we have 
shown that the presence of strongly polarizable co-solvent 
molecules in the diffuse layer (at the large values of surface 
potential) leads to a more long-ranged electrostatic poten¬ 
tial profile. This effect results from the fact that with an 
increase of the local co-solvent concentration in the diffuse 
layer the local dielectric permittivity increases. The latter 
leads to a decrease of the electrode charge screening. 

In conclusion, we would like to discuss the weak point of 
the present theory and its possible generalization. First, 
our model should be valid at the low co-solvent concen¬ 
tration only. Second, within the present study we have 
neglected the excluded volume interactions for the ions 
as well as for the co-solvent molecules. Nevertheless, as 
Kornyshev clearly showed in the work 19 , the excluded 


volume interactions should be also important in the re¬ 
gion of large surface potential for a correct evaluation of 
the differential capacitance. Moreover, it is evident in ad¬ 
vance the account of the excluded volume interactions for 
the co-solvent should decrease the discussed effect of the 
increase of the differential capacitance. However, we be¬ 
lieve that the latter is still significant even taking into 
account the excluded volume of the co-solvent molecules 
as well as the electrolyte ions. A discussion of the impor¬ 
tance of the excluded volume contributions together with 
the polarizability and permanent dipole of species on the 
properties of the electric double layer is a subject of the 
future publications. 
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